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Abstract 

The Landauer-Butikker formalism is an important formalism to study mesoscopic systems. Its 
validity for linear transport is well established theoretically as well as experimentally. Akkermans 
et al [Phys. Rev. Lett. 66, 76 (1991)] had shown that the formalism can be extended to study 
thermodynamic properties like persistent currents. It was earlier verified for simple one dimensional 
systems. We study this formula very carefully and conclude that it requires reinterpretation in quasi 
one dimension. This is essentially because of the presence of evanescent modes in quasi one dimension. 

PACS numbers: 







I. INTRODUCTION 



Due to the technological advances in nano-fabrication, it is possible to realize such small 
systems that the quantum mechanical coherence length of the electron extends through out 
the length of the sample. Quantum interference phenomena strongly determines the thermo- 
dynamic and transport properties of these so called mesoscopic systems^. 

Mesoscopic phenomenon can occur in canonical systems as well as in grand canonical sys- 
tems. Mesoscopic systems are so small that even measuring probes (like voltage probe and 
current probe) can make the system a grand canonical systent^. A canonical mesoscopic sys- 
tem is well described by the Hamiltonian of the isolated system but that is not the case for grand 
canonical mesoscopic systems. A grand canonical mesoscopic system, by definition, is coupled 
to a reservoir with which it can exchange electrons. Mesoscopic systems are so small that the 
reservoir can drastically change the states of the system and this has to be explicitly accounted 
for-. One can take this into account by solving the Schrodinger equation of the leads and 
the system as a scattering problem. This approach is essentially known as Landauer-Butikker 
formalism. This formalism is thus different from the way we deal with large grand-canonical 
systems with the help of the grand partition function. 

For example, if an Aharonov Bohm flux is applied through the center of a ring, the ring gets 
magnetized and a persistent current is generated in the ring^. This current arises because of 
the vector potential that changes the phase of the wave function in the ring and is another form 
of Aharonov Bohm effect which is an interference phenomenon, whereas the magnetization is 
a thermodynamic property is a consequence of that. One can cite many similar phenomenon^. 
Persistent currents has been studied for more than a decade, theoretically as well as exper- 
imentally. If the ring is isolated then the persistent current is carried in some eigen-states. 
Whereas if it is open and connected to reservoirs, then persistent current is carried in reso- 
nant and non-resonant states that are typical of scattering states. Persistent current has been 
studied in several such grand canonical systems like a ring connected to a single reservoir or 
to many reservoirs^. If it is connected to many reservoirs at different chemical potentials, then 
non-equilibrium currents can co-exist with equilibrium persistent currents. In such open sys- 
tems several interesting effects have been predicted, like current magnification in the presence 
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of transport^'^'^, directional dependence of persistent currents-, current magnification effect in 
equilibrium systems in absence of transport current^, etc. 

In order to realize a mesoscopic grand canonical system we connect the ring to reservoirs 
that are at fixed chemical potentials as is schematically shown in Fig. [TJ The left reservoir has 
a chemical potential /ii and the right one has a chemical potential fi2- The reservoirs can also 
be at a finite temperature T. The ring is threaded by an Aharonov-Bohm flux. If fii > fi2, then 
there is a transport current (which is a non-equilibrium current) through the regions I and II. 
There is however no transport current in the ring. The ring will carry a persistent current which 
is an equilibrium current. Thus in the present geometry, the equilibrium persistent currents 
and non-equilibrium transport currents are spatially separated. 

The Landauer-Buttiker approach proposes that an equilibrium phenomenon like persistent 
current in such an open system as that in Fig. 1, can be obtained from solving the scattering 
problem^. Akkermans et al^ related the persistent current IS to the S matrix by the following 
formula. 

1 c^log[det(g)] 

Such a simple mathematical relation between the persistent current inside the ring and the S- 
matrix obtained from the wave-function far away from the ring is rather novel and resulted in a 
flurry of theoretical activities^. While it is established (from theoretical and experimental point 
of view) that conductance (a non-equilibrium phenomenon) can be obtained from the S-matrix, 
Akkermans approach may prove to be the first step to obtain any equilibrium phenomenon from 
the S-matrix, that is a step towards a mesoscopic version of fluctuation dissipation theorent^. 
The correctness of Eq. [T]has been explicitely verified in one dimension (ID) but not in quasi 
one dimension (QID). That complexities arise in QID due to the presence of evanescent modes 
has been observed recently^>^ although it is not very well known to the community^. For 
example Friedel sum rule^ and Buttiker-Thomas-Pertre formulaM breaks down in the presence 
of evanescent modes. So in this work we undertake the task of verifying if Akkermans formula 
is valid in QID. In case of Akkerman's formula, in this paper, we can show analytically how the 
evanescent modes complicates things. Earlier works^^ii^ on different formulas are essentially 
numerical verifications. Persistent current in the geometry of Fig. 1 has been studied earlier. 
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but always using the wave function. The S-matrix was never used and comparison was not made 
between Akkerman's approach and the usual wave function approach. This is done explicitely 
for the first time in QID in this work to show that the ID result cannot be extended to QID 
due to the presence of evanescent modes. 



II. MODEL AND METHOD 



As shown in Fig. 1, we consider a ring coupled to a wire. The scattering solution for this 
geometry is discussed in detail in our earlier work^. Here we outline some points with respect to 
calculating the RHS of Eq. [T]which was not done earlier. There is a (5-potential impurity present 
in the ring at any arbitrary position X [Fig[T]. We apply Aharonov-Bohm flux (p through the 
ring, perpendicularly to the plane of the paper. We consider two modes of propagation because 
it will show the shortcomings of Akkerman's formula and that can be generalized analytically 
to any number of modes. The Schrodinger equation for a QID wire in presence of a 5-potential 
at X = 0, y = Hi is (the third degree of freedom, i.e. 2;-direction, is usually frozen by creating a 
strong quantization^) 

' ^^)^Vc{v)\'^[x,y) = E^[x,y) (2) 



2m* dx^ d'lp- ' 

where the x-coordinate is along the wire and the y-coordinate is perpendicular to the wire. 
Here m* is the electron mass and E is the electron energy. The wave-function in a ring can be 
obtained by solving the above equation with periodic boundary condition where we assume the 
ring to be so large that its curvature can be neglected. Here Vc{y) is the confinement potential 
making up the quantum wires in Figure [H The magnetic field just appears as a phase of \E'(x, y) 
that will be accounted for while applying boundary conditions. Eqn. [2] can be separated as 

~ 2m* ~ 2m* ^ ^ 

and 

+ V,{y)]xn{y) = E^xn{y) (4) 



2m* dy'^ 

Here we consider that electron is propagating along x direction. This means in regions / and // 
of Fig. 1, X direction is along the arrows. In region J J J, the x direction is along the line joining 



P and Q. And in regions IV and V, the x direction is along the perimeter of the ring. One 
can choose different axes in the different regions as the matrix equations for mode matching is 
independent of this choice^. The confinement potential Vc{y) in different regions is then in the 
y (transverse) direction. It can be seen from Eqs. [H] and H] that 



E = En + 



2m* 



(5) 



We take Vc{y) to be a square well potential of width W, that gives Xn{y) = sin[^(2/ + ^)]. So 
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2^^^. Hence, in the first mode, 



2m*E vr^ 



is the propagating wave-vector and in the second mode 
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2m*E 47r2 



(7) 



is the propagating wave-vector. We have chosen 2m* = 1 and h = 1. 

When electrons are incident along region I (in FigH]) in the first mode the scattering problem 
can be solved exactly. The solution to Eqn. [3] in region / becomes 
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Similarly, in region //, ///, IV and V we get 
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where r[-^, r[2, g'u and g[2, A, B, C, D, E, F, G, H, J, K, L and M are to be determined by 
mode matching. 

Note that at P and Q we have a three legged junction that is schematically shown in Fig. 
2. In a previous work^ we proposed a form of junction scattering matrix Sj for a two channel 
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FIG. 1: A ring connected to an infinite wire. A 6 function potential is present in the ring at position X. 
A chemical potential difference (/ii — ^^2) between the left reservoir (LR) and the right reservoir (RR) 
drives a transport current through the regions I and II. The ring is pierced by an Aharonov-Bohm 
flux that drives an equilibrium current called persistent current in the ring. 



quantum wire that can be easily generalized to any number of channels. For the 6 potential 
impurity at X we use the scattering matrix S}, that was derived by Bagwell^. One can match 
the wavefunctions and conserve the currents by using these S'-matrices that give us a set of 
linear euations. We calculate the coefficients A, B, C, D, E, F, G, H, J, K, L and M 
numerically by matrix inversion. 

Persistent current can be computed from the wave-function. 




FIG. 2: A 3-leg junction. 



w 




2 



(13) 



5 



Here the index ki implies that this is the current due to an incident electron in ki channel on 
the left. Similarly currents are generated due to incident electron in ki channel on the right, 
/c2 channel on the left and k2 channel on the right. So the net observable persistent current is 



IW = 2JPy('=i) + 2JPy('=2) 



(14) 



From Eq. [T3]we get 



_ 2/o(|Ep - + |Gp - \H\ 



2^(fcl) 



(15) 



where Jq = 



We also calculate the scattering matrix elements r[^, r[2,g'ii, g'12: ''^21; ^'22? fi'21 ^"7 matrix 
inversion to form the scattering matrix of the system as 



S 



(16) 



^ r'li g'u 9'i2 

^21 ^22 921 922 

9ii 9i2 '"ii ^'12 

\ 921 9-22 ^21 ^"22 / 

Substituting (TH] in [H Eq. [T] too can be written as a sum of four termsi^, where each term 
consist of scattering matrix elements due to incidence in a particular momentum channel. That 
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where 
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Although not imphed by the notation, the currents defined above (Eqs. [TH and [T7j) are 
actually differential currents in an infinitesimal energy range dE. The integration of these will 
give the actual measurable currents. The integration limits depend on the chemical potential fii 
and /i2. Temperature can be included through Fermi function. All expressions so far is derived 
for both modes being propagating. Earlier it was shown that IW = IS for one dimensional ring 
coupled to a reservoir-ii. We shall show below that when all modes are propagating then one 
gets IS = IW, but not when we include evanescent modes. This is because when evanescent 
modes are present then some expressions can be analytically continued to include evanescent 
modes but not all of them. 

III. INCLUSION OF EVANESCENT MODES 

E is the energy of incidence that can be varied as an external parameter by tuning the 
chemical potentials of the reservoirs. When < 2m*EW^/h^ < An'^, then it can be seen 
from Eq. [7] and Eqs. 8-12 that k2 mode becomes evanescent. Eqs. 8-12 are still solutions to 
Schrodinger Eq. [S] implying electrons in the ring can be coupled to an evanescent channel due 
to scatteringi^. A single impurity can couple an electron to the evanescent second channel. 
Scattering at the junctions can also couple an electron to the evanescent second channel. No 
electron can be incident from (or emitted to) oo along the evanescent 2nd channel. So the 
scattering problem has to be solved with an incident electron in ki mode in the left lead and 
an outgoing electron in ki mode in the right lead. Hence the S matrix becomes 2x2 and is 
given by 

\ g'n r'n 

Although the S matrix is 2 x 2, its calculation has to be done by using the 6x6 junction S-matrix 
Sj that is defined in Ref. p] and the 4x4 S-matrix 5*^ for the 6 function potential that is defined 
in Ref. p^]. This is essential because evanescent modes can be coupled inside the ring without 
violating any physical principle like conservation of energy. g[2, ^'12 ^^^^ non-zero, but 

they do not carry any current. They are not S-matrix elements any more. Rather they define 
the coupling to evanescent modes. Unitarity should imply | ^ + | p= 1 and indeed 
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FIG. 3: IW/Iq and IS/h vs Sir'^m* EW"^ /h'^ . The system parameters are ^ = I2 = k = 1, yi = 0.1, 
a = P = 0.3 and 7 = 4. 



we get this from the junction matrix defined by Sj. This imphes that Sj is appropriate for 
accounting for reahstic multichannel situations. 



IV. RESULTS AND DISCUSSIONS 



In a real system, there are always propagating modes as well as evanescent modes. As is 
evident from Eqs. E] and [TJ evanescent modes have higher transverse energy than the propa- 
gating modes. Also higher the n value of the evanescent mode, higher is its transverse energy. 
This energy is at the cost of the propagation energy which becomes more and more negative 
for higher n evanescent modes. There will be a natural cut off as very high energies cannot be 
realized in a quantum wire. In our simplistic approach we will first consider a case when there 
are two modes in the wire, both being propagating. We will then consider a case when one mode 
is propagating and the other is evanescent. In the first case we will verify that Akkerman's 
approach gives exactly the same result as the current calculated from the wave-function (i.e., 
IS = IW). In the second case there are complexities. It will be argued that such complexities 
will persist in a real system where there will be many evanescent modes. 

When both modes are propagating: First we consider the energy range 
9tt^ (i.e. 39.478 < 2m*EW^/h'^ < 88.826). Substituting this E in Eqs. Eland [7] we can see that 
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FIG. 4: IS^^^^/Iq and IW'^'^^^/Iq versus 87r2m*£;w7/i^- The system parameters are h = I2 = I3 = 1, 
y. = 0.1, a = /3 = 0.3, 7 = 4. 

both the modes are propagating. The nature of current is plotted in Fig [31 The figure shows 
that the current IS obtained from Akkerman's formula (that is from ^-matrix) is identical 
with IW (that is from wave function). In fact, IS^'^'^^ and IS^'^^^ are individually identical with 
JVr^'^i) and IW^^^\ respectively. These are shown in Fig. H] and Fig. [5l This implies that 
JiyC^i) is the same algebraic expression as IS^^'l Similarly for IW^''^^ and IS^''^\ And also IS 
and IW are the same algebraic expressions. 

When one mode is evanescent: Now consider the energy range vr^ < 2m*EW'^/h'^ < 47r^ (i.e. 



9.87 < 2m*EW^/h^ < 39.477) so that ^2 = \ - becomes imaginary (^2 ^ ^2) while 



evanescent mode. Evanescent mode current can be calculated by directly applying Eqn. [T3] to 
evanescent mode wave functions or it can be calculated by analytically continuing propagating 
mode current [Til to below the barrier (that is ^2 ^^2)- We have already argued that IS and 
IW are the same algebraic expression. Under the transformation k2 iK,2, applied to both IS 
and IW, they definitely remain the same algebraic expression. However, Akkerman's formula 
takes a different meaning in this regime where there are evanescent modes. This is essential 





In this regime the ring contains one propagating and one 
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FIG. 5: IS^^^^/Io and /VF^'^s)//^ versus 87r'^m*EW'^ /h"^ . The system parameters are h = I2 = I3 = 1, 
y. = 0.1, a = /3 = 0.3, 7 = 4. 

because the Akkerman's formula in Eq. [1] is related to the S-matrix and the transformed 
expression {ISk2~>iK2) cannot be obtained from the S-matrix. We know that no electron can 
be incident along the evanescent channel. So IW^'^'^^ and IS^''^^ are zero. So IW is now just 
equal to IW^'''^\ where k2 has been analytically continued. If one assumes that in this regime 
IS is equal to JS'^'^^^ where k2 is analytically continued, then obviously once again IS = IW, as 
they are the same algebraic expression. However, one can see that JS''-'^^^ (see Eq. [T8l) cannot 
be obtained by substituting the S-matrix (Eq. [20]) . into Akkerman's formula (Eq. [1]). If we do 
this substitution, then we will get the first term and the third term in Eq. [THl (where of course 
^2 — > iK,2 transformation is taken care of). We will not get the 2nd and 4th terms as 
g'12 are not S-matrix elements any more. In fact, we see from Fig [6] that the difference between 
IS and IW is quite large. 

Since IS is now reduced to just two terms, one may ask the question that will it give 
the partial current in the propagating channel only. Because this partial current also consists 
of two terms only. Note from Eqn. [IS] that the total current JiyC'i) = IW^''^^ + IW^''^\ 
where IWi'"'^ = 2/o(| E p - | F p)^^!), E and F being the wave function amplitudes in the 
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FIG. 6: IS/Iq (solid line) and IW/Iq (dashed line) versus Svr^ m*E'Ty^//i^ of the system when second 
channel is evanescent. The system parameters are li = I2 = h = !> Hi = 0.1, a = (5 = 0.3, 7 = 4. Here 
IS/Iq is obtained by substituting 5 given in Eq. [50] into Eq. 1. 
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FIG. 7: IS/Iq (solid line) and IW^ /Iq (dashed line) versus ^Ti'^m* EW"^ /h? . The system parameters 
are h = h = h = !> Hi = 0.1, a = (3 = 0.3, 7 = —3.7. Here IS/Iq is obtained by substituting S given 
in Eq. [50] into Eq. 1. 

propagating channel and IW2^^^ = 2Iq{\ G {"^ — \ H p)('^'i), G and H being the wave function 
ampUtudes in the evanescent channel We have plotted IWi'^^^ in Fig. [3 The figure shows that 
IS diff'ers from IWi''^\ So it is confirmed that IS obtained from Eq. [T] neither give the total 
measurable current of the system nor the partial current through the propagating channel. 
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V. CONCLUSIONS 



For realistic mesoscopic rings connected to leads, there are always evanescent modes. The S- 
matrix is always defined by the propagating modes only. For such systems one cannot directly 
apply Akkerman's formula. Instead one should start with a model where all the modes are 
made propagating. One should apply Akkerman's formula to the S-matrix of this system and 
then analytically continue this expression for the current to the situation where an appropriate 
number of modes are evanescent. While the Landauer-Buttiker approach is still inevitable as 
the evanescent modes are obtained due to an incident electron that is scattered to evanescent 
modes, the formula given in Eq. 1 is no longer strictly valid in presence of evanescent modes. 
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